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Abstract 

We investigate the source imaging for a granular pion-emitting source model in high energy 
heavy ion collisions. The two-pion source functions of the granular sources exhibit a two-tiered 
structure. Using a parametrized formula of granular two-pion source function, we examine the two- 
tiered structure of the source functions for the imaging data of Au-|-Au collisions at Alternating 
Gradient Synchrotron (AGS) and Relativistic Heavy fon Collider (RHIC). We find that the imaging 
technique introduced by Brown and Danielewicz is suitable for probing the granular structure of the 
sources. Our data- fitting results indicate that there is not visible granularity for the sources at AGS 
energies. However, the data for the RHIC collisions with the selections of 40 < centrality < 90% 
and 0.20 < kq^ < 0.36 GeV/c are better described by the model with granular emission than that 
of one Gaussian. The model with granular source has more parameters than the simple Gaussian, 
hence can describe more complicated shapes. 
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I. INTRODUCTION 



Two-pion Hanbury-Brown-Twiss (HBT) interferometry is a powerful tool of detection 
of the space-time structure of particle-emitting sources produced in high energy heavy ion 
collisions [l, 2, 3, 4|. In conventional two-pion HBT analysis one needs fitting the two-pion 
HBT correlation functions with parametrized formulas to obtain quantitatively the source 
space-time results. So the explanations of the HBT results are model depended. Imaging 
technique introduced by Brown and Danielewicz ^, y, Q] is a model-independent way. It can 
be used to extract the source geometry pictures (source function) directly from the two-pion 
correlation functions. This technique has been developed and used in high energy heavy ion 



collisions 
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i2,ii3i,y,ii5i, 
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19|,l20|. 



Recently, there has been much progress in understanding of the process of nucleus-nucleus 
collisions at RHIC. However, there are still many unsolved problems. One of them is the 



so-called HBT puzzle, i?out /-Rside 
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24| . Here i?out and -RsiHp are the transverse 
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26|. In Ref. |27| 



HBT radii parallel and perpendicular to the pion pair momentum 
a granular source model of quark-gluon plasma (QGP) droplets was put forth to explain 
the HBT puzzle. The suggestion was based on the observation that in the hydrodynamic 
calculations for the granular source the average particle emission time scales with the initial 
radius of the droplet, whereas the spacial size of the source is the scale of the distribution 
of the droplets. For a granular source with many of the small droplets distributed in a 
relatively large region, the HBT radius i?out can be close to -Rside In Ref. 28] the 

authors further investigated the elliptic fiow and HBT radii as a function of pion transverse 
momentum for an improved granular source model of QGP droplets. They argued that 
although a granular structure was suggested earlier as the signature of a first-order phase 



transition 
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37l | , th e occurrence of granular structure may 
not be limited to first-order phase transition 28|, l38|] . The large fluctuations of initial matter 

42] in high energy heavy ion collisions may facilitate the occurrence 



distribution 
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4G, 
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of instability of the system during its subsequent expansion and fragmentation to many 



28 



granular droplets together with surface tension effect 

Recent researches on event-by-event two-pion Bose-Einstein correlations in smoothed 
particle hydrodynamics indicate that the particle-emitting sources produced at RHIC en- 
ergy are inhomogeneous and there is a granular structure of many "lumps" (droplets) j43|. 
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For a simple granular source model we will show that the two-pion source function has a 
two-tiered structure. In small relative coordinate r region, the source function exhibits an 
enhancement because of the higher density in the droplets. Previous RHIC experimental 
imaging researches are mainly focused on the long-range tail of the source functions at large 
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18, M, uSl, and the deviations of the source function from Gaussian distribution in 



the large r region are believed mainly the contribution of long-lived resonances |17l . |44| |. In 
this paper we will focus our attention on the source functions in small r region. We will 
investigate the imaging of granular sources. We will examine the two-tiered structure of 
the source functions for the AGS and RHIC imaging data of Au+Au [lO, 113]. Our results 
indicate that the imaging technique is suitable for probing the granular structure of the 
particle-emitting sources. There is not visible granularity for the sources in Au+Au col- 
lisions at AGS energies. However, the data for the RHIC collisions with the selections of 
40 < centrality < 90% and 0.20 < k-r < 0.36 GeV/c are better described by the model with 
granular emission than from that of one Gaussian. 



II. IMAGING TECHNIQUE 



For the convenience of discussion later in the paper^we start out with a brief review of 
the imaging technique of Brown and Danielewicz 5|, |6|, |7|. 



Based on the Koonin-Pratt formulism 
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461, the two-pion HBT correlation function 



may be expressed in the center-of-mass frame of the particle pair as 



15|: 



C(q)-1 = JdvKici,r)S{r), 



(1) 



where q = pi — p2 is the relative momentum of the pion pair, r is the relative separation 
of emission points of the two particles, K{c{,r) = |$q(r)p — 1, where $q(r) is the relative 
wave function of the pair. Neglecting the final-state interaction of the pion pair, one has 



«q-r/2\ 



(2) 



In Eq. ([T]), >SYr) is the so-called two-particle source function. It may be written with Wigner 



function as 



S{r) = jdt jd^RdTD(R+r/2,T+t/2,pi) 
xD(R-r/2,T-t/2,p2). 



(3) 



Here the Wigner functions are normalized particle emission rates, 

For a spherically symmetric source function, ^(r) = S{r), performing the angle integra- 
tions on the right hand of Eq. ([1]), one gets the angle-averaged version of Eq. ([1]) as, 

7^(ginv) = C(ginv) - 1 = 47r y drr^Kiqi^,, r)S{r). (5) 



Here ginv = - qh 

K{qinv, r) = sin(ginvr)/ (gmv?") . (6) 

Equation now is suitable in any frame and the problem of imaging becomes inverting 
K{qinv,r) with measured correlation function 7^(ginv)- After expanding the source function 
S{r) in 6 — spline basis [3], 

S{r) = J2S,B,{r), (7) 
i 

the inversion problem reduces to the problem of solving the matrix equation 

niq.) = J2K^,S,, (8) 

j 

where 7l{qi) denotes the value of the correlation function at the ith bin of ginv, and 

'qi-Aq/2 



Kij = — rfginv / drr'^K{qi^^,r)Bj{r), (9) 

^1 Jqi-Aq/2 Jo 



where Aq is the bin size of ginv 

In present paper, the minimization package MINUIT 0] was used to minimize the 
between the measured and calculated correlation functions. 



III. SOURCE FUNCTION IMAGING FOR STATIC GRANULAR SOURCES 



In this section we examine the source function imaging for static granular sources. Al- 
though static source is not a realistic case, its source function may still be used as a 
parametrized formula in the analysis for evolving sources. 
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For the granular source model, particles are emitted from dispersed droplets 27, 



28 



35l . I36l . |37[. Assuming that the granular source has the same droplets and the 



distribution of the particle emission points in a droplet has Gaussian form, the normalized 
source density distribution is given by [361] 

N , 



DM 
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A^(V27ra)3 



1=1 



2a2 



(10) 



where a is the "radius" of the droplets and Xj is the spatial coordinate of the ith droplet 
center. Inserting this distribution into Eq. ([3]), we obtain 



S{r) 



iV2(V47ra)3 
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X exp 
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V 4a2 J 


I 2a2 ; 



(11) 



where X 



X,; 



Xj, a is the angle between r and Xj,-. The source function presents 



fluctuation due to the factor exp 



r|Xij I cos a 
2^ 



Due to limited number of produced particles per event, conventional two-pion HBT in- 
terferometry are analyses based on averages over events. The correlation functions are 
obtained from the correlated pion pairs (the two identical pions in each of the pairs are from 
the same event) of all sample events. Accordingly, the source function for the mixed-events 
of the granular source, S'^''^"(r), should be the average of Eq. (ITT]) over all the events. 
Assuming that the droplet centers in the granular source obey the Gaussian distribution, 
P(X,)~exp(-Xf/2i?2j,weget 
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Granu / 
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: 6xp 
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4(a2 + Rl 



(12) 



Compared 

B,[io,Q, 



Compared with the two-particle source function of the Gaussian source model usually used 

1 



S 



Gauss / 



47ri?ga)3 



(13) 



the source function of the granular source consists of two exponential terms, which corre- 
spond to that the two particles from the same droplet and from different droplets respec- 
tively. This "two-tiered structure" of the source function is consistent with the two-tiered 
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FIG. 1: The two-particle source functions of granular sources. 



structure of the correlation function for the granular source [30|, 13^. It is most apparent 
when N = 2, and disappears when oo. For R^^. ^ a^, the coefficient ratio of the two 

terms in Eq. is {R^/af/{N - 1), which is a characteristic quantity for the two-tiered 
structure. Considering also > 2 for the granular source model, we introduce the quantity, 



(14) 



to characterize the granularity of the sources. For granular sources the values of ^ are 
between (0,oo), and a source will have not granularity when ^ < 0. In Fig. 1 we show 
the curves of S'~^'^^^^{r) as a function of the droplet number for the granular sources with 
i?gr = 4.5 fm and a = 1.5 fm. One can see that the the two-tiered structure is obvious for 
finite droplet numbers of the granular sources. 

We next examine the imaging of the static granular source. In Fig. 2(a) we show the 
two-pion correlation function (CF) obtained from 2 x 10^ simulated two-pion events for the 
granular source with Rgj- = 4.5 fm and a = 1.5 fm. For comparison, the correlation function 
for a Gaussian source with i?ga = 4.5 fm is presented in Fig. 2(b). In the simulations, pions 
are emitted thermally from the sources at freeze-out temperature Tf = 150 MeV. Figure 
2(a') and (b') show the two-pion source functions (circle symbols) extracted from the two- 
pion correlation functions by the imaging technique. The curves in Fig. 2 (a') and (b') are 
the results of the source function fit (SFF) with the formulas XS '^'^'^'"^^{r) and XS '^'^^^^{r). 
Here A is the parameter of source coherent factor in HBT interferometry 

mm 

and 
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FIG. 2: (a, b) The two-pion correlation functions for granular and Gaussian sources, (a', b') The 
two-pion source functions extracted by imaging technique for the granular and Gaussian sources. 

^Granu^^-j and 5* Gauss giveii by Eqs. (|T2l) and f|T3l) . respectively. The corresponding 
fitting results are also presented in the figure. For the granular source the x^/NDF for 
the granular SFF is 0.22, which is much smaller than that of 3.88 for the Gaussian SFF. 
The curves in Fig. 2(a) and (b) are the restored correlation functions (RCF) calculated by 
Eq. ([5]). From Fig. 2 one can see that although the two-pion correlation functions for the 
granular and Gaussian sources are almost the same in shape, the granular source function 
exhibits a clear enhancement in small r region, which reflects the higher source density in a 
droplet. So imaging technique is suitable for probing the granular structure of the sources. 

Figure 3(a), (b), (c), and (d) show further the source functions for the granular sources 
with various source parameters. It can be seen that the two-tiered structure is more obvious 
for smaller droplet number and smaller droplet radius a. One can still observe the two- 
tiered structure even ^ ~ 2. 



IV. TWO-TIERED STRUCTURE OF SOURCE FUNCTIONS FOR EVOLVING 
SOURCES 

In this section we investigate the two-tiered structure of the source functions for evolving 
sources. First we consider a simple granular source model of evolving QGP droplets. We 
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FIG. 3: The source functions for the granular sources with various parameters. 



assume that all of the droplets in the source have the same initial conditions and evolve 



hydrodynamically in the same way 



[27 



STLAu 



suggested by QCD lattice gauge results [48|, |49|, l50|| with the transition temperature Tc = 165 



equation of state of the entropy density 



MeV and the temperature width of the transition AT = O.OSTc is used in the hydrodynamical 
calculations 27, In our calculations the initial energy density of the droplets is 



taken to be eo = 3.75 TcSc, which is about two times of the density of quark matter at Tc 
0, [5^. The initial distribution of the droplet centers is given by a Gaussian distribution 
with standard deviation Rq. For the case with an additional collective radial expansion, the 
droplet centers are assumed to have a constant radial velocity in the center-of-mass frame 

37|. The source freeze-out temperature is taken to be Ty = 150 



27 



of the granular source 
MeV. 

Figure 4(a) and (b) show the source functions S{r) (circle symbols) obtained by the 
imaging technique from the two-pion correlation functions C(ginv) for the granular sources 
with the parameters Rq = 5.0 fm, Vd = 0.5, and A = 1. The droplet number and initial 
droplet radius for the granular source of Fig. 4(a) are 5 and 2.5 fm, and they are 15 and 1.5 
fm for the granular source of Fig. 4(b). One can see that the source functions of the granular 
sources have enhancements in small r region. With the granular SFF results one gets that the 
values of ^ for the two granular sources are 4.85±2.94 and 4.66±2.83, respectively. The large 
^ values indicate large granularity for the sources, which is consistent with the observations. 
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FIG. 4: The source functions for the evolving granular sources [(a) and (b)] and the expanding 
Gaussian sources [(c) and (d)]. 

For comparison, in Fig. 4(c) and (d) we show the source functions obtained by the imaging 
technique from the two-pion correlation functions for the sources with Gaussian distribution 
and additional radial expanding velocities 0.3 and 0.6, respectively. The standard deviation 
for the Gaussian distribution is taken to be 5.0 fm. From Fig. 4(c) it can be seen that for 
the Gaussian distribution source with smaller expanding velocity the granular and Gaussian 
SFF curves are almost overlapped. The granular SFF gives a very large A^. For a large 
one can see that the first term of the fitting formula Eq. (fT2|) approaches zero and the 
two-tiered structure disappears. One can also see from Eq. ( |T2i) that in this case the fitting 
formula is almost independent of and the parameters Rg,- and a can be reduced to one 



parameter ■\/R^^ + cfl as i?ga in Eq. ( |T3l) . So the fit is insensitive to the parameters 
and -Rgi-. From Fig. 4(d) one can see that in small r region the two SFF curves are almost 
overlapped and there is only a small difference between the two SFF curves in large r region. 
We find that the granular SFF result of A^ is less than 2. In this case the source has not 
granularity. With the granular SFF results for the two expanding Gaussian sources we get 
that the values of the characteristic quantity of granularity ^ are 4.813 x 10~^ ± 0.002 and 
— 12.900 ± 4.185. They indicate that there is not granularity for the sources. 

We next examine the two-tiered structure of the two-pion source functions for the imaging 



data in the Au+Au collisions at AGS 



and RHIC 



17|. The circle symbols in Fig. 5(a) 



(b), (c), and (d) show the two-pion source functions for 2, 4, 6, and 8 AGeV Au+Au 
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FIG. 5: The two-pion source functions (o symbols) for 2, 4, 6, and 8 AGeV Au+Au collisions 
and the SFF curves. 



collisions, respectively lOj. The circle symbols in Fig. 6(a) and (b) show the source functions 
for ^/sNN = 200 GeV Au+Au collisions with different cut conditions of centrality and 
average transverse momentum /ct of the pion pair. In Figs. 5 and 6, the solid and dashed 
curves are our granular and Gaussian SFF curves. The fitting results are listed in Table I. 
It can be seen that the xVNDF for the granular SFF are very small. This indicates that 
the data errors are large for the granular SFF and the granular SFF may distinguish more 
complicated source shapes if there is enough statistics. From Fig. 5 one cannot observe the 
two-tiered structure of the source functions in small r region. When we use the granular 
parametrized formula XS'^'^'^'°^^{r) to fit the source functions we find that the errors of are 
the same order of the values of A^. In Table I we present the fit results for fixed A^ = 50. 
One can see that the values of ^ in Table I for the collisions at AGS energies are very small. 
So there is not visible granularity for the sources. 

From Fig. 6 it can be seen that for the higher /ct most central collisions (a) the source 
function has not obvious two-tiered structure and the spheroidal (dot-line) 17|], granular, 
and Gaussian SFF curves in small r region are almost overlapped. It indicates that there 
is not visible granularity for the source. However, for the lower peripheral collisions (b) 
one can find an obvious two-tiered structure of the source function. In small r region the 
source function has an enhancement relative to the Gaussian SFF curve fsee the insert in 
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FIG. 6: The two-pion source functions (o symbols) for y^J/vTV = 200 GeV Au+Au collisions 
and the SFF curves. 



TABLE I: The results of granular and Gaussian source function fits (SFF). 



AGS(a) AGS(b) AGS(c) AGS(d) RHIC(a) RHIC(b) 



Granular 


i?gr (fm) 


6.31 ±0.23 


5.92 ±0.18 


SFF 


a(fm) 


2.67 ±0.16 


2.92 ±0.18 




A 


0.93 ± 0.05 


0.77 ±0.02 




N 


50 (fixed) 


50 (fixed) 




X'/NDF 


0.04 


0.05 






0.28 ±0.08 


0.17 ±0.05 


Gaussian 


i?ga(fm) 


6.26 ±0.11 


6.22 ±0.07 


SFF 


A 


0.84 ±0.03 


0.73 ± 0.02 




X'/NDF 


0.62 


0.46 



5.28 ±0.19 


5.39 ±0.23 


5.23 ±0.21 


4.53 ±0.10 


2.51 ±0.21 


2.35 ±0.17 


3.90 ±0.13 


2.56 ±0.05 


0.57 ±0.02 


0.65 ±0.04 


0.24 ±0.01 


0.39 ±0.01 


50 (fixed) 


50 (fixed) 


3.54 ±0.61 


4.61 ±0.35 


0.01 


0.06 


0.35 


0.71 


0.19 ±0.07 


0.25 ± 0.09 


1.56 ±0.97 


2.12 ±0.55 


5.53 ± 0.09 


5.41 ±0.12 


5.05 ±0.13 


3.79 ± 0.03 


0.53 ±0.02 


0.59 ± 0.03 


0.21 ±0.01 


0.32 ±0.01 


0.32 


0.41 


3.72 


10.87 



Fig. 6(b)]. Based on the granular source explanation, the enhancement in small r region 
indicates that there are small droplets with higher density in the particle-emitting source. In 
Table I, the value of ^ for the case (b) is larger than that for the case (a), which is consistent 
with the observations. Further investigation for the reasons of the enhancement of source 
function in small r region will be of great interest. 
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V. SUMMARY AND CONCLUSION 



We investigated the source imaging for a granular pion-emitting source model. The two- 
pion source functions of the granular sources exhibit a two-tiered structure, which can be 
characterized by the quantity ^ = {Rgr:/a)^/{N — 2). In small relative coordinate r region, 
the granular two-pion source functions have an enhancement because of the higher density 
in the droplets. We find that the imaging technique is suitable for probing the granularity 
of the pion-emitting sources. Using a parametrized formula of granular source function, we 
examine the two-tiered structure of the source functions for the imaging data of Au+Au 



collisions at 2, 4, 6, 8 AGeV [10|] and ^/sNN = 200 GeV [l7l]. Our analysis results indicate 
that there is not visible granularity for the sources produced in the collisions at the AGS 
energies and at RHIC energy with the selections < centrality < 20% and 0.48 < < 0.60 
GeV/c. However, the data for the RHIC collisions with the selections 40 < centrality < 90% 
and 0.20 < < 0.36 GeV/c are better described by the model with granular emission than 
from that of one Gaussian. The model with granular source has more parameters than the 
simple gaussian, hence can describe more complicated shapes. 

Although our granular parametrized formula of source function is obtained from a static 
granular source model and does not include the effect of droplet overlap, the fitting results 
with the formula for evolving sources have still referential meaning. In this paper we only 
examine one-dimension imaging of granular sources. Because the longitudinal dynamics 
at the RHIC energy is very different from that at AGS energies, the examinations of the 



source imaging in different directions 



3, 



15 



16 



44| and at intermediate energies {i. e. 



SPS energies) will be of great interest. Further investigations on the source granularity and 
its variation with the centrality and particle transverse momentum in collisions will be also 
interesting issues. 
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